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RESUMEN 

El analisis de una estructura aislada requiere de un modelo matematico que 
proporcione curvas fuerza-desplazamiento que permitan idealizar el comportamiento 
no lineal del aislador. La propuesta tiene como finalidad usar un modelo simple 
utilizando las propiedades del mismo y validarlo con datos experimentales de un 
proyecto real, que permita contar con un modelo que estima en forma razonable el 
comportamiento del aislador Se presenta un modelo de analisis sismico de 
aisladores de triple pendulo de friccion, el mismo que es validado con el ensayo de 
81 aisladores, sometidos cada uno a tres ciclos de carga, con una duration de 12 
segundos. Por otra parte, se analiza la variation de: coeficiente de rozamiento, 
rigidez efectiva, factor de amortiguamiento y perfodo de vibration, encontrados en 
forma experimental y con el modelo propuesto en los 81 FPT. 

Palabras Claves: Aisladores de Triple Pendulo de Friccion. 

ABSTRACT 

The analysis of an isolated structure required a mathematical model to 
provide lateral force - displacement curves allowing to idealize the nonlinear 
behavior of the device. The proposal aims to use a simple model using its 
mechanical properties and validate it with experimental data of a real project, 
allowing to have a model that estimates in a reasonable manner the behavior of 
insulator. We present a model for seismic analysis of a triple friction pendulum, 
which is validated with the trial of 81 devices, subjected each to three cycles, with a 
duration of 12 seconds. On the other hand, is analyzed the variation of coefficient 
of friction, effective stiffness, damping factor and vibration period, found 
experimentally and with the model proposed in the 81 FPT. 
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INTRODUCCION 


Terremotos recientes han demostrado que, incluso modernos codigos han 
limitado los danos en elementos estructurales, pero la mayor concentration de 
perdidas se concentra en los danos en los componentes y elementos no 
estructurales, asf como en la pedida de funcionalidad (Zayas 2013). Ante esta 
realidad, es importante considerar sistemas como el aislamiento de base que 
permite limitar el dano estructural y en componentes y elementos no estructurales, 
permitiendo con esto que las estructuras logren niveles de desempeno optimos 
reduciendo las perdidas. (Aguiar etal. 2008, Kawamura et al. 2000) 

Los dispositivos de aislamiento de base, los cuales principalmente se han 
concentrado en dos tipos como son los llamados elastomeros y friccionales 
(Naeim and Kelly, 1999). Los primeros en este caso, los elastomericos han sido 
desarrollados e implementados inicialmente existiendo principalmente tres tipos de 
aisladores elastomeros y son los conocidos como: Low Damping Rubber Bearing; 
High Damping Rubber Bearing y Lead Rubber. (Constantinou 2012). 

De igual manera los friccionales se pueden clasificar en tres tipos: Simple 
Concave Friction Pendulum Bearing; Double Concave Friction Pendulum Bearing; 
Triple Concave Friction Pendulum Bearing. La investigacion continua y es asf 
como acaba de aparecer los denominados: Quinto Pendulo de Friccion. (Lee y 
Constantinou 2016) 

Cabe mencionar que pese a las ventajas y a las aplicaciones existentes 
(Chistopupoulus 2006), existen limitaciones en su aplicacion o limites, 
principalmente para estructuras muy esbeltas, con una gran cantidad de pisos, por 
el efecto P - A. Ademas de lo relacionado a la componente vertical que afecta a 
elementos no estructurales como son los cielos rasos; este tema ha sido 
investigado en el laboratorio en Japon, E-Defense (2011). 

Nos centraremos en los aisladores friccionales de triple pendulo de friccion 
FPT, ya que seran los que se coloquen en los nuevos Centros de Investigacion, de 
la Universidad de Fuerzas Armadas-ESPE. (UFA-ESPE). Estos aisladores 
combinan superficies de friccion con fuerzas de restauracion por la geometrfa de 
los mismos (Fenz 2006). 

Los aisladores friccionales dobles y triples llamados de segunda y tercera 
generacion respectivamente, tienen algunas ventajas, por ejemplo son mas 
compactos, capaces de adaptar su desempeno en relation a la demanda, 
aumentar la capacidad de desplazamiento, reducir la velocidad en el movimiento, 
con lo cual se evita variation en los coeficientes de friccion. Otro aspecto notable 
del aislador de segunda generacion con respecto al de primera es que se reduce 
las respuestas estructurales mejorando el desempeno de los componentes y 
elementos no estructures. (Fenz y Constantinou 2008). 

Los aisladores de triple pendulo de friccion estan conformado por un 
aislador interior de radios de curvatura R 2 ,R 3 , dentro de otro aislador cuyas placas 
exteriores tienen radios R V R 4 . De tal manera que se tiene dos aisladores en uno 
solo. Esto permite que tengan menores dimensiones con respecto a los de primera 
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y segunda generation y que tengan una mayor capacidad de desplazamiento. 
(Constantinou et at. 2016). 

En la UFA-ESPE se estan construyendo 6 Bloques Estructurales con 
aisladores FPT8833/1 2-1 2/8-6, el mismo que se indica en la figura 1. En total se 
necesitaron 81 aisladores los mismos que fueron ensayados considerando tres 
cargas verticales de 250, 200 y 1 00 T. (EPS 2015). 



Figura 1 Geometrfa del aislador FPT8833/1 2-1 2/8-6, utilizado en la Universidad de 
Fuerzas Armadas ESPE, de Ecuador. 

En la figura 2 se muestra la Implantation General del Centro de 
Investigation y Posgrado de la Universidad de las Fuerzas Armadas (ESPE) con el 
HATCH se encuentran los edificios con aislacion sismica y sin HATCH los edificios 
con disipacion. 



ESC HA 


Figura 2 Implantation General del Centro de Investigation y Posgrado de la 
Universidad de las Fuerzas Armadas (ESPE). 
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En la figura 3 se muestra el avance general de cada bloque estructural, el 
proyecto esta en su fase de montaje de estructura metalica y construccion de 
accesos principals la obra presenta un avance general del 40 %. 



c) 


d) 



Figura 3 a) Estado Actual del Bloque 1 montaje de estructura metalica b) Estado 
Actual de Bloque 2 y 3 montaje de aisladores en primera planta c) Estado Actual 
del Bloque 4 y 5 construccion de accesos d) Estado Actual del Bloque 6 y7 
montaje de correas e) Estado Actual del bloque 8 colocacion de losa tipo deck. 
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MODELO DE TRES FASES 


Los radios de curvatura de las placas exteriores e interiores del aislador 
FPT pueden ser diferentes al igual que las alturas h t con lo que la capacidad de 
desplazamiento d t son diferentes (para i = 1:4); con lo que se tiene 12 
condiciones geometricas y cuatro condiciones de race fa diferentes en cada una 
de las placas. En este caso el modelo de cinco fases propuesto por Fenz y 
Constantinou (2007, 2008); Fadi y Constantinou (2010) es el mas adecuado. 

Ahora para el caso del aislador FPT8833/1 2-1 2/8-6. Las condiciones 
geometricas se reducen a 6 debido a que los radios de curvatura de las placas 
exteriores son iguales y lo propio sucede con los radios de las placas interiores; 
algo similar se tiene con las alturas y la capacidad de desplazamiento como se 
observa en la figura 1. Se puede tener dos coeficientes de race uno para las 
placas exteriores y uno para las placas interiores. Para estas condiciones Me Vitty 
y Constantinou (2015) propusieron un modelo de tres fases para definir las curvas 
de histeresis que relacionan el desplazamiento horizontal con la fuerza de code. 
Las ecuaciones que definen este modelo, son: 


Ri.eff - Ri~ hi Para i de 1 a 4. 


( 1 ) 


dl = 


Ri.eff 

Ri 


( 2 ) 


Donde R t es el radio de curvatura; hi es la altura, presentada en la figura 
1 ; R ief f es el radio de curvatura efectivo; d* es la capacidad de desplazamiento. El 
subindice i, varfa de 1 a 4. A continuation se indican las 3 fases o regimes del 
modelo. 


• Regimen I 

El desplazamiento se da solamente en las placas 2 y 3. 

0 < u < u* 

(3) 

(4) 

Donde u es el desplazamiento lateral del aislador; F es la fuerza lateral 
aplicada; W es el peso que gravita sobre el aislador. A la izquierda de la figura 4 
se aprecia al aislador interior desplazandose en las superficies 2 y 3; a la derecha 
se muestra el diagrama de histeresis correspondiente. 


u* = 2 (jUi - jU 2 ) R 2i e ff 

W 


F = 


2 R 


2 ,eff 


U+ /J. 2 W 
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Regimen 1 


Displacement 

Figura 4 Desempeno de aislador en Regimen I. 
Fuente: Me Vitty y Constantinou (2015) 


• Regimen II 

El movimiento en el aislador interior Mega a los topes y empieza a 
desplazarse en las superficies 1 y 4. Normalmente es en este regimen en que 
trabaja un aislador ante un sismo de intensidad moderada y alta. Las ecuaciones 
se indican a continuacion y en la figura 5 se presenta la curva de histeresis 
correspondiente. 


u* <u < u** 
xT = u* + 2 dl 
W , 

F = — (u-u*)+ li t W 

Z K 1 ,eff 


(5) 

( 6 ) 



Displacement 


Figura 5 Desempeno de aislador en Regimen II. 
Fuente: Me Vitty y Constantinou (2015) 


• Regimen III 

Este regimen se presenta cuando el sismo es extremadamente fuerte y el 
aislador interior se topa con los topes exteriores, en estas condiciones empieza a 
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deslizarse el aislador interior en las superficies 2 y 3. Las ecuaciones de calculo y 
el desempeno se indican en la figura 6. 

u — n — 

u cap = 2d{ + 2d* 2 


( 7 ) 

( 8 ) 



Displacement 


Figura 6 Desempeno de aislador en Regimen III. 

Fuente: Me Vitty y Constantinou (2015) 

MODELO PROPUESTO 

El modelo que se propone es exclusivo para la fase II, donde va a trabajar 
el aislador, por lo regular. Pero puede ser aplicado tambien para la fase I. Cabe 
indicar que difiere del modelo indicado por propuesto por Me Vitty y Constantinou 
(2015), primero que no se considera un valor de resistencia a cero 
desplazamiento igual a Qd que es funcion de la friccion y del peso y tambien en 
la descarga no se considera la llnea vertical de longitud 2 Q d , en lugar de ello se 
considera una llnea inclinada como se vera posteriormente. 

En la figura 7 a) se presenta el modelo propuesto por Me Vitty y 
Constantinou (2015) para la fase II, se aprecia que cuando se inicia la descarga lo 
hace con una llnea vertical indicada en el parrafo anterior y luego con una llnea 
inclinada cuya pendiente es la misma que la rigidez elastica. 

Ahora, se propone en la figura 7 b) que cuando se inicia la descarga se lo 
haga directamente con una rigidez igual a la rigidez elastica. Es decir que se 
trabaje con un modelo bilineal cuyo coeficiente de friccion viene definido por la 
siguiente ecuacion. 


= Ai - (i“i - Mz) 


Rlef_ 

Rlef 


0 ) 


Donde |Ui,/r 2 , son los coeficientes de roce de las placas exterior e interior; 
es el coeficiente de roce equivalente. 
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Displacement Displacement 

a) b) 

Figura 7 Modelos para la fase dos: a) propuesta de Me Vitty y Constantinou (201 5) 
y b) modelo propuesto en este estudio. 

Las ecuaciones que definen modelo bilineal, son las siguientes: 


<T _ ^ H' 

(10) 

1 w 


fcr, = 

(11) 

F = gW + k p q 

(12) 

k ef = ~ 

(13) 

rj, o / W 

T = 2 n 

(14) 

Aj 9 k ef 


Donde W es la carga vertical que gravita en el aislador; q es el 
desplazamiento lateral en el aislador que se lo obtiene en forma iteractiva; f e(? es 
el factor de amortiguamiento equivalente; k e f es la rigidez efectiva o rigidez 
secante; T es el perfodo del aislador; g es la aceleracion de la gravedad. 

RESULT ADOS EXPERIMENTALES 


En la figura 8 a) se observa parte de los 81 aisladores FPT adquiridos por 
la UFA-ESPE a EPS, en la 8 b) el transporte de 4 de ellos en un monta carga 
hacia el area del ensayo; en la 8 c) se observa un aislador, sin la protection 
exterior durante el ensayo y en la 8 d) la curva de histeresis que relaciona los 
desplazamientos con las fuerza lateral en los ciclos tres ciclos de carga que 
tuvieron una duration de 12 segundos y un desplazamiento lateral que se halla 
alrededor de las 12 pulgadas. 
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c) d) 

Figura 8 a) Parte de los isladores de UFA-ESPE; b) transporte de los aisladores 
para ser ensayados; c) ensayo de un aislador; d) relaciones desplazamiento fuerza 

y ciclos de carga aplicados. 

Los aisladores por efecto del movimiento (figura 8 b), ocasionan que la 
base rfgida y el aislador interior se muevan, de tal manera que no se encuentran 
centrados para lograr que esten completamente en la posicion original, 
manualmente se asigna un ciclo de carga lo que se obteniendo la grafica de la 
figura 9 a), con Ifneas que luego son eliminadas, lo propio sucede con la parte final 
del ensayo hasta que vuelva a su posicion inicial con lo que el desplazamiento 
lateral es cero, figura 9 b). 


Non lineal Behavior Isolator 40 



Displacement (cm) 


a) 



b) 
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Horizonlal Displacement (mm) 

c) 

Figura 9 a) Curva inicial que contempla cojinete descentrado; b) Curva en la que 
se ha eliminado la carga dada manualmente; c) aproximacion de un modelo 
numerico a los resultados experimentales. 

Relation Cortante expresado en funcion de la carga vertical vs desplazamiento 

horizontal de aislador 40. 

Finalmente, se encuentra una curva que mejor se ajuste a los tres ciclos 
de carga, en base a la cual se determina los coeficientes de roce, siguiendo el 
modelo de cinco fases de Fenz y Constantinou (2007, 2008). Para el aislador 
FPT8833/1 2-1 2/8-6 en la figura 8 se presenta la curva de histeresis 
correspondiente. 

En la figura 10, f t corresponde al coeficiente de roce del aislador interior 
que se ha denominado /^ 2 ; / 2 ,/ 3 , a los coeficientes de roce de las placas exteriores 
denominadas Experimentalmente se determinan estos coeficientes de roce 
y de acuerdo al modelo de cinco fases, EPS (2015), determino la rigidez efectiva, 
k ef , el factor de amortiguamiento equivalente % eq y el perfodo de vibration T, 
asociados al desplazamiento lateral de 12”. 
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SHEAR 



Figura 10 Curva de histeresis code desplazamiento de aislador FPT8833/12- 

12/8-6 

Fuente: EPS (2015) 


RESULTADOS 

Con los coeficientes de roce, obtenidos experimentalmente, EPS (2015) 
encontro la rigidez efectiva, el amortiguamiento equivalente y el perfodo del 
aislador asociados al desplazamiento de 12”, utilizando el modelo de 5 fases. 
Ahora en este estudio se encuentran los parametros indicados, empleando el 
modelo bilineal propuesto. 

EPS (201 5) reporta los coeficientes de roce en cada ciclo de histeresis y el 
valor promedio, para cada uno de los 81 aisladores sfsmicos FPT. Es importante 
indicar que 61 aisladores fueron ensayados con una carga vertical de 250 T. A 
otros 10 aisladores se los ensayo con una carga vertical de 200 T. A los 10 
restantes se les aplico una carga de 100 T. 

Se obtuvieron tres tipos de curvas, una para valores promedios del 
coeficiente de roce, otra para valores maximos y otra para valores mmimos. 
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5.1 Coeficientes de Friccion. 

En la figura 1 1 se compara los valores de Coeficientes de Friccion cuando 
la carga vertical es de 250 T. Los valores hallados experimentalmente son 
producto de los 61 ensayos realizados por EPS. 


FRICTION COEFFICIENTS OF OUTER PLATES FOR 250 Ton 



a) 


FRICTION COEFFICIENTS OF INNER PLATES FOR 250 Ton. 



b) 


Figura 11 Comparacion de coeficientes de friccion para 250 T de carga vertical, a) 
Coeficientes de friccion de las placas exteriores ( %): valores medios, valores altos 
y valores minimos b) Coeficientes de friccion de las placas interiores ( u 2 ): valores 
medios, valores altos y valores minimos 

En la figura 12 se compara los valores de Coeficientes de Friccion cuando 
la carga vertical es de 200 T y 100 T. Los valores hallados experimentalmente son 
producto de los 20 ensayos realizados por EPS. 



a) 



b) 


Figura 12 Comparacion de coeficientes de friccion para 200 T y 100 T de carga 
vertical, a) Coeficientes de friccion de las placas exteriores ( u x ): valores medios, 
valores altos y valores minimos b) Coeficientes de friccion de las placas interiores 
(u 2 ): valores medios, valores altos y valores minimos 
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5.2 Rigidez efectiva 

En la figura 13 se compara la rigidez efectiva cuando la carga vertical es 
de 250 T hallada con valores: maximos, mmimos y promedios de coeficiente de 
friccion. Se encontro que los valores hallados experimentalmente son ligeramente 
mayores a los encontrados con el modelo bilineal propuesto. La mayor diferencia 
con los dos modelos es menor al 4%. 



Figura 13 Comparacion de rigidez efectiva, para 250 T., de carga vertical, a) con 
valores promedios de race; b) con valores maximos de coeficiente de race y c) con 
valores mmimos de coeficiente de race. 


En las figuras 14 se comparan la rigidez efectiva, cuando la carga vertical 
es de 200 T. (graficos de la izquierda) y 100 T. (graficos de la derecha), 
respectivamente y se encuentra un comportamiento similar, aunque para algunos 
casos la rigidez efectiva que se halla es igual con los dos modelos. 
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Figura 14 Comparacion de rigidez efectiva, para 200 T. (izquierda) y 100 T. 
(derecha), de carga vertical, a) con valores promedios de race; b) con valores 
maximos de coeficiente de race y c) con valores mmimos de coeficiente de race. 

5.3 Factor de amortiguamiento equivalente 

En las figuras 15 y 16 se presentan los factores de amortiguamiento 
hallados para carga vertical de 250 T. (61 ensayos); 200 T. (10 ensayos) y 100 T. 
(10 ensayos). Se observa que el factor de amortiguamiento equivalente 
encontrado con el modelo propuesto es ligeramente mayor que el encontrado en 
forma experimental. 
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Figura 15 Comparacion de factor de amortiguamiento equivalente, para carga 
vertical de 250 T., de carga vertical, a) con valores promedios de roce; b) con 
valores maximos de coeficiente de roce y c) con valores minimos de coeficiente de 

roce. 


Cato a) W=200 T. Average Friction Values 
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Figura 16 Comparacion de factor de amortiguamiento equivalente, para carga 
vertical de 200 T. (izquierda) y 100 T. (derecha), de carga vertical, a) con valores 
promedios de roce; b) con valores maximos de coeficiente de roce y c) con valores 

minimos de coeficiente de roce. 
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5.3 Periodo de vibracion 

En la figura 17 se presentan los perfodos del aislador cuando la carga 
vertical es de 250 T., y en la figura 16 cuando la carga vertical es de 200 T. 
(izquierda) y de 100 T. (derecha). En el siguiente apartado se comenta sobre los 
resultados obtenidos. 


'Case a) W=250T. Average Fri ction Value 


■■ EqE ii i e il U l 


UhhMlL 


lM* 
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Caso b) W=250T. Maximum Friction Values. 
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Figura 17 Comparacion de factor de amortiguamiento equivalente, para carga 
vertical de 250 T., de carga vertical, a) con valores promedios de race; b) con 
valores maximos de coeficiente de race y c) con valores mfnimos de coeficiente de 

race. 
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Figura 18 Comparacion de Perfodos, para carga vertical de: 200 T. (izquierda) y 
100 T. (derecha), de carga vertical, a) con valores promedios de roce; b) con 
valores maximos de coeficiente de roce y c) con valores mmimos de coeficiente de 

roce. 


VARIACION DE LOS RESULTADOS OBTENIDOS 


Dos puntos interesan comparar con los resultados obtenidos, el primero la 
variacion que tiene la rigidez efectiva, factor de amortiguamiento equivalente y 
perlodo encontrados experimentalmente y segundo comparar esos resultados con 
los que se hallan con el modelo propuesto. Para este objetivo en las tablas 1 a 3 
se presentan los valores medios y desviacion estandar con los datos de las figuras 
13 a 18. 
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Tabla 1 Variacion de la rigidez efectiva encontrada con valores: Promedio, 
Maximos y Mmimos para tres condiciones de carga vertical. 


Valores 

Modelo 

W = 250 T. 

W = 200 T. 

W = 100 T. 



Kf 

a Kef 

Kf 

° K ef 

Kf 

° K ef 



T 

T 

T 

W 

T 

T 

W 

T 

Promedio 

Experimental 

101.49 

3.12 

83.02 

1.15 

45.91 

0.97 


Propuesto 

99.92 

2.96 

81.67 

1.02 

44.62 

1.13 

Maximos 

Experimental 

106.32 

3.11 

86.88 

1.37 

47.61 

0.92 


Propuesto 

105.51 

3.49 

85.87 

1.07 

46.57 

0.76 

Mmimos 

Experimental 

97.95 

2.92 

79.98 

1.22 

44.69 

1.02 


Propuesto 

95.60 

3.26 

78.45 

1.50 

43.18 

1.56 


Tabla 2 Variacion del factor de amortiguamiento encontrado con valores: 
Promedio, Maximos y Mmimos para tres condiciones de carga vertical. 


Valores 

Modelo 

w = 

250 T. 

W = 200 T. 

w = 

100 T. 



? 


? 


? 


Promedio 

Experimental 

0.25 

0.0103 

0.2634 

0.0043 

0.29 

0.0062 


Propuesto 

0.27 

0.0109 

0.2864 

0.0045 

0.31 

0.0089 

Maximos 

Experimental 

0.27 

0.0089 

0.2788 

0.0036 

0.30 

0.0042 


Propuesto 

0.29 

0.0116 

0.3037 

0.0044 

0.33 

0.0051 

Mmimos 

Experimental 

0.24 

0.0099 

0.2527 

0.0053 

0.28 

0.0093 


Propuesto 

0.26 

0.0135 

0.2718 

0.0073 

0.30 

0.0133 
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Tabla 3 Variacion del Periodo encontrado con valores: Promedio, Maximos y 
Minimos para tres condiciones de carga vertical. 


Valores 

Modelo 

w = 

250 T. 

w = 

200 T. 

w = 

100 T. 



T 

(5.) 

G j* 

is.) 

T 

is-) 

G j 

is.) 

T 

is.) 

G j 

is.) 

Promedios 

Experimental 

3.12 

0.0443 

3.08 

0.0222 

2.93 

0.0297 


Propuesto 

3.17 

0.0471 

3.14 

0.0195 

3.00 

0.0384 

Maximos 

Experimental 

3.17 

0.0499 

3.14 

0.0220 

2.97 

0.0349 


Propuesto 

3.08 

0.0509 

3.06 

0.0191 

2.94 

0.0239 

Mmimos 

Experimental 

3.04 

0.0441 

3.01 

0.0226 

2.88 

0.0465 


Propuesto 

3.24 

0.0559 

3.20 

0.0306 

3.05 

0.0554 


COMENTARIOS Y CONCLUSIONES 

Cabe indicar que el modelo bilineal propuesto es consistente y proporciona una 
estimacion de la respuesta de la estructura, que podria ser compatible con los 
comentarios de Me Vitty y Constantinou. 

El modelo propuesto es validado con los datos experimentales proporcionados 
por EPS, en base al aislador friccional utilizado en el Nuevo Centro de 
Investigaciones de la Universidad de Fuerzas Armadas, Espe. 

Los valores de Rigideces, Amortiguamiento y periodos del modelo experimental y 
propuesto con valores maximos, minimos y promedio permite observar que el 
modelo analitico bilineal es compatible con el experimental. 
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